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Abstract 



Recent experiments on deformation of semiconductors show an abrupt 
the variation in the critical resolved shear stress r Y with tempera ure 
?SESL a change in the deformation mechanism at a cn«cal tonpe^ure 
T In the cases examined so far in our laboratories (Case Western Reserve 
University and Poitiers) and elsewhere, this critical temperature appears to 
coincide approximately with the brittle-to-ductile transition ^P^^^ 
™ this paper, new deformation experiments performed on the wide-bandgap 
semtondJrfor 4H-S1C over a range of temperatures at two stram rates are 
bribed together with a transmission electron microscopy characterization of 
S I disbcations below and above T c . Based on these, and results recendy 
reported on a few III-V compound semiconductors a new model for the 
dTformation of tetrahedrally coordinated materials at low and high 
tem^ratures is. proposed, and the relation of the transition m deformadon 
mode to the transition in fracture mode (brittle to ductile) is discussed. 

§1. Introduction 

The plasticity and fracture of semiconductors have interested materials scientists 
for a long time because, compared with possibly any other material, semiconductors 
can be grown as crystals of exceptionally high quality. In particular, since the inven- 
tion of the transistor in 1949, techniques have been developed to grow many semi- 
conducting crystals to contain a very low density of dislocations (or none in the case 
of silicon) or, for that matter, any other type of defect. Consequently, cracks or 
dislocations can be introduced in an initially defect-free (or nearly so) material in- a 
controlled manner and their behaviour studied under the application of a tensile or 
shear stress. * 



1 1 . Plastic properties of semiconductors 
Some of the early work on the plasticity of semiconductors was carried out by 
Alexander, Haasen and their co-workers and was reviewed in their seminal paper 
(Alexander and Haasen 1968). More recent reviews of the work in this area are by 

§ Email; pxp7@cwru.edu 
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The topics reviewed in these papers <^ «* J ^ ot defor- 

nis^te—^d^ 
of the high Peierls potential in «h crysttk, the ^ ^« " ™ 

are dissociated according to the following reactions: 

|{110) -i(121)+i(2U), (1) 

$<2lTo)-i (lioo) +|(ioio), 

S5 Sen £2 l^^e^ -d Suzuki 1999). In addition in the caseof 

1 comp^carbon core). We shall later argue that the core nature of J-drffcn* 
parZ * such dislocation loops plays a significant role m the mode of deformauon 
that occurs in tetrahedrally coordinated materials. 

13 Dynamics of dislocation glide 
In general, the dislocation velocity v in a tetrahedraUy coordinated material can 
be expressed in the following form: 

/ &Gl{r t T)\ (2 ) 
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i(g) 1^=4- [1010] 
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Figure 1. A dissociated hexagonal loop 011 the (0001) slip plane in SC. 

where the stress r and temperature T dependence of the dislocation velocity are 
Tncorporated in the free-energy barrier AG(r,r) fo V^%^ 0 V£ 
absence of applied stress (r = 0), the free-energy barrier A<?„ = AG(0, T) is 
assumed to have the same temperature dependence as the shear modulus M . In 
many' cases, over a large range of stresses and temperatures, the dislocation velocity 
can be expressed in terms of the following phenomenological equation: 



v=x (£f exp <r 



kT 



(3) 



where AJ? B is a stress- and temperature-independent activation enthalpy change for 
dislocation glide, and A, r 0 and the stress exponent m are constants (A includes an 
activation entropy term for dislocation glide and t 0 = 1 MPa if t is m megapascals). 
In this form, the stress dependence of the dislocation velocity appears in the pre- 
exponent and the activation enthalpy is assumed to be a constant A# g . 

The relation between macroscopic and microscopic plastic properties of a crystal 
are given by Orowan's equation for the plastic strain rate 7: 



j = pbv, 



(4) 



where p is the density of mobile dislocations and b is the magnitude of the Burgers 
vector of dislocations producing the shear strain 7. Thus, combining equations (.4) 
and (2), one has 



7 = pbvo exp 



^ AG T (r,r) ^ 



(5) 



In this case, the stress dependence of AG T = A// r - T A5 T (where AH T and AS r are 
the enthalpy and entropy changes, respectively) can be determined from the activa- 
tion volume, V* = -(dG T /(h) T = -(dHr/dr) T , using stress relaxation tests or 
strain-rate jump tests (for example Kubin (1974), Kocks et al. (1975)). 
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If equation (3) is applicable in the range of stresses and temperatures under 
consideration, then one can determine the ««^^t ^ 
Aft and the stress exponent m from a direct measurement of dis location velocity 
^function of ten^rature and applied stress. A different way to 
activation parameters of dislocation glide involves measurement of the JJ^S? 
as a function of temperature and applied stress. This is again obtamed by combining 
the macroscopic parameters of yield and microscopic properties of d^location 
through substituting equation (3) for dislocation velocity in Orowan s relat.on (4) 
to give (Alexander and Haasen 1968) 



r Y = 7 1 /" exp 



(6) 



Thus it is possible to determine the constant Q by measuring the critical resolved 
shear stress as a function of temperature at a fixed strain rate. In addition, the 
constant n can be determined by performing strain-rate jump tests (for example 
Samant and Pirouz (1998)). The constants n and Q are related to the phenomeno- 
logical parameters^ and Ai/ g describing the dislocation velocity m equation (3). 
The relation between the constants n and and AH S and Q can be derived from the 
kink diffusion model of dislocation glide as n = m + 2 and A# g - nQ (for example 
Alexander and Haasen (1968)). 

The facts that dislocations in semiconductors are straight when produced at low 
temperatures and he along the <110> directions are highly suggestive that their 
motion is controlled by a high Peierls stress and that, at finite temperatures, they 
move by the mechanism of kink pair formation and kink migration. Employing the 
kink diffusion model of Hirth and Lothe (1968), the activation enthalpy Atf g can be 
related to the enthalpy Affe required for kink pair formation and the enthalpy 
Atf km required for kink migration. Although, formerly, it was thought that, just 
as in fee crystals, the kink migration energy is rather negligible and Ai/ g « Af^r, 
this assumption has recently been questioned and the kink migration energy in 
tetrahedrally coordinated materials may turn out to be large and even 
Aif^ > Affjtf . For details of measuring A# g , Ai/^ and Atf^, reference can be 
made to the aforementioned review papers of George and Rabier (1987) and Rabier 
and George (1987). 

1. 4. Fracture behaviour of semiconductors 
The work on the fracture of semiconductors started much later than the work on 
the plasticity of these materials. At low temperatures, semiconductors are brittle 
materials; the elemental semiconductors (silicon, germanium and diamond) cleave 
preferentially on a {111} plane whereas the (cubic) compound semiconductors pre- 
ferentially cleave on a {1 10} plane (similarly the hexagonal semiconductors cleave on 
a (1010} plane). One of the earliest papers on fracture of semiconductors reported 
measurement of the surface energy of {111} planes in silicon by cleavage (Gilman 
1960). However, the first systematic fracture study of a semiconductor was made by 
St John (1975) who performed careful fracture measurements at different strain rates 
on pre-cracked silicon single crystals over a wide range of temperatures (between 
-196 and 1000°C). He found a very narrow temperature range (a few degrees 
Celsius) over which silicon changed from a brittle to a ductile behaviour. This 
brittle-to-ductile transition (BDT) temperature T B dt w *s found- to be a sensitive 
function of the applied strain rate e (or, equivalent^, the rate K of application of 
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the stress intensity factor K) and was thermally activated with an activation energy 
close to that for dislocation glide. The seminal work of St John on silicon was later 
reoeated and extended by a number of researchers including extensive studies by 
Brede and Haasen (1988), Michot (1988), Haasen et al (1989), Hirsch et al 
(1989a b) Hirsch and Roberts (1991), Hsiah and Argon (1994), George (1998). 
These researchers extended the work of St John to include the influence of various 
variables on the BDT temperature, such as electrically inactive impurities (e.g. oxy- 
gen and hydrogen), electrically.active dopants (e.g. donors and acceptors) and crack 
olane orientation. Tests were also conducted on a few other semiconductors such as 
germanium (Serbena and Roberts 1994), GaAs (Michot et al 1988, Michot and 
George 1989) and SiC (Maeda and Fujita 1989). Most of the experiments mentioned 
have been performed at a constant rate K of increase in the stress intensity factor K 
under a mode I loading geometry- In this mode, the tensile stress a is applied 
perpendicular to the crack plane (see, for example figure 10 later). The microscopic 
studies of the fractured samples, by X-ray topography (Michot and George 1982, 
1986, George and Michot 1993) or by transmission electron microscopy (TEM) 
(Chiao and Clarke 1989), have shown that there is hardly any dislocation 
activity below J BDT , while an avalanche of- dislocations nucleate at this critical 

temperature. • . . 

The most quantitative model proposed for the BDT in silicon is by Hirsch 
and Roberts (1991) where the shielding of the crack front by dislocations emitted 
from there competes with the rise in the stress intensity factor K to the critical 
value K Xc . The dislocations are assumed to be emitted from a few special sites 
along the crack tip; these primary ' sites are separated from each other by an 
average distance d x . Subsequently, the loops emitted from these sites help to 
activate nucleation from other, more difficult sites along the crack front. These 
secondary sites are separated from the primary sites by an average distance d 2 
(usually taken to be equal to 4/2). The main feature of the model is that it is 
only when the. emitted dislocations shield every point of the crack front that the 
material becomes ductile and the crack front is not able to propagate catastro- 
phically. Under a constant with an increasing K t the shielding is considered to 
be effective if the dislocation loops expand at a sufficiently rapid rate that they 
cover the whole of the crack front before K reaches X lc . In this sense, the 
mobility of dislocations plays the major role in the Hirsch-Roberts model because 
the loops from the primary sources must traverse a distance d { - d 2 to cover 
every point on the crack before K increases to the critical value K lQ . The dis- 
tances d{ and d 2i as well as the stress, intensity factors and K 2 , required to 
nucleate dislocation loops at the primary and secondary sites respectively are used 
as fitting parameters in the model. With the help of these four fitting parameters, 
reasonably good agreement with experimental results have been obtained for the 
prediction of J B dt as wel1 as for , ^ variation in T Bm with K\ for more detail, 
refer to the original papers of Hirsch and Roberts (1991, 1997). 

The Hirsch-Roberts model has recently been criticized by Khanta (1994) and 
* Khanta et al. (1994a,b) on the grounds that a thermally activated process cannot 
explain the sharpness of T^ot (see, however, the series of correspondence between 
Hirsch and Roberts (1996a,b) and Khanta and Vitek (1996)). In analogy with the 
Kosterlitz-Thouless (1973) model of paraelectric-to-ferroelectric transition, Khanta 
et al have proposed a model based on autocatalytic nucleation of dislocations in the 
material at T BDT (see also Sun et al (1998)). As will be seen below, the sharpness of 
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the BDT temperature in our model arises from the nucleation of an avalanche of 
^pSdislocations under the applied stress aided by thermal acUvafon. 

1 5 Plasticity of semiconductors in the brittle regime 
Many interesting experiments on the yield behaviour of semiconductors involve 
defection at low Temperatures where the crystals are usually brittle. A number of 
suSTxperiments have b'een performed by Vickers Indentation where 
stress field has a hydrostatic component that prevents the fracture of the material 
Superimposed on this hydrostatic component, there are large tensile or C0 ™P^e 
and shear stresses that plastically deform the material by shear (for example, Hirscn 
et al (1985)). However, it is not usually possible to obtain quantitative data on me 
yield stress and dislocation glide activation parameters from indentation tests iirus 
the more useful low-temperature deformation experiments involve uniaxial tests sucn 
as tension or compression. One way to perform such tests is to pre-oeform tne 
material at a high temperature (in the ductile regime) and subsequently to deform 
it (usually by compression) at a lower temperature (in the brittle regime) (for exam- 
ple Wessel and Alexander (1977)). It should be mentioned that, if the initial disloca- 
tion^density in the material is sufficiently high, a pre-deformation is not essential and, 
as long as the deformation is performed in a carefully aligned compression jig under 
very low stram rates, it is possible to deform the sample plastically at temperatures 
down to a few hundred degrees Celsius below the BDT. A different technique for low- 
temperature deformation is to perform the compression (or even tensile) test m the 
presence of a superimposed hydrostatic pressure in order to prevent the fracture of 
the test sample. An early study in this respect was that of Castaing et al. (1981) who 
deformed silicon in the temperature range 300-600°C in a Griggs apparatus, usually 
used for deformation of brittle geological materials. Rabier et al. (1985), Demenet 
(1987) and Demenet et al. (1989) later performed similar types of experiments on 
silicon and on GaAs. More recently, Suzuki et al. (1998, 1999a,b) and Edagawa et al. 
(2000) used a similar technique to deform InP, InSb, GaAs, and GaP; as before, the 
hydrostatic pressure prevented fracture of the deformation samples before they plas- 
tically yielded. . . 

In our laboratories (Case Western Reserve University (CWRU) and Laboratoire 
de Metallurgique Physique (LMP)), research has concentrated on two different hex- 
agonal polytypes of the wide-bandgap semiconductor SiC As part of his PhD thesis, 
Samant (1999) investigated the plastic behaviour of 4H-SiC and 6H-SiC by compres- 
sion tests over a wide range of temperatures and strain rates. It should be mentioned 
that, because single crystals of these two semiconductors were not available until 
recently, prior to the work of Samant, there were only a few reports of deformation 
tests on single-crystal 6H-SiC and almost none on 4H-SiC. Some important studies 
that existed were compression temperature range by Fujita et al. (1987) on Acheson- 
grown 6H-SiC crystals over the temperature range 1 300-1 600°C, and creep tests on 
Cree-grown 6H-SiC by Corman (1992). Samant's experiments covered the tempera- 
ture range 550-1400°C at strain rates varying from 6.5 x 10~ 4 to 3.1 x 10" s . He 
was able to deform the crystals at temperatures below their usual range of BDT 
because of the rather high initial density of dislocations in the samples (about 
10 3 -10 4 crn~ 2 ) ) careful alignment of the sample in the deformation jig and with the 
use of very low strain rates (Samanz and Pirouz 1998, Samant et al. 1998, Samant 
1999). More recently, Demenet et al. (2000) repeated Samant's experiments on 4H- 
SiC (grown at Cree Research, Inc.) obtaining much more data in the temperature 
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orvY_n*n=C at both CWRU and LMP. Here we report some of the experi- 
range 900-1360 C at . d discuSS ^ implications for the plas- 

mental measurements of Demenet et at. anu ^ 
ticity and fracture of semiconductors. 

§ 2. Yield behaviour of 4H-SiC single crystals 

2.1 Experimental procedures 
The details of the experimental procedure for d eformingttie *C 
been given elsewhere (Samant 1999). Brief* a ^S^S^X^ 
at Cree Research , Inc was 0™^^^^^^ faction. From the 

aces The^mpression tests were made over the temperature . range ; 800-13*0 C 
v£g two strain rates varying by an order of magnitude <£ « 3-6 x 10 and 
T^urVh The series of deformation tests at a strain rate ol 
illdxlO-V 1 Were performed at CWRU, while those at a stram rate of 
- 2 6 x 0-V 1 were performed at 1PM; no hydrostatic confining pressure was 
Ufa W 0* bese " 1S P tt should also be mentioned that tests m the J*»Habom. 
useoinany »J 10 -5 s -i basically reproduced the results obtained earuei at 
CWRU tfor"mor de.ails . « X (2000)). Following .he mechanical 
« s to fo^raUd to the (0001) slip planes were prepared from the samples 
Ttl'Jn^J, ^Lratures above and below the transition temperature. In addition 
S^rSSen^o characerizethe But^-^onhed^caUons 
the laree-anele convergent-beam electron diffraction (LACBED) technique ^nerab 
"stS luSToSu e, al. 1992) was used to identify .he nature of drslocauon 
cores (Ning and Pirouz 1996). 



2.2. Results 3 



Figure 2 shows a plot of In (r Y ) versus l/T for two strain rates: e « ^ 1 u 
and 2 6 x 10- 6 s" 1 . As in the previous experiments on 4H- and 6H-&C (Samant and 
1998 Samant et al. 1998, Samant 1999), there is a sharp transUton at a cri^a 
temperature T c . A significant point about these carves is the shift of T, to aJowe 
emoerature as the strain rate is lowered Thus, whereas T C «1100C tor 
?~?6 xTo'V 1 the transition temperature is shifted by 70<C to a lower tempera- 
u7e o I « IWC when the slower strain rate e « 2.6 x l£ £» uj«Ut shouid 
be noted that while macroscopic plastic deformation of 4H-SiC (and 6H-S1C) is 
^lan' ely easy at T> T„ it becomes very limited at T < T< and, when the sample. 
£ ^ not Med catastrophically, deformation is very often accompa me b h 
appearance of .extensive microcracking. These observauons from a fafly p a 
crystal with no evidence of any cracking when deformed at T > T t to a crystal 
S undergoes very limited plasticity accompanied by profuse 
when deformed at T < T c , indicate that the BDT ^%^ r f^mi) and 
close to T c . This is also consistent with the experiments of Fujita et al. (1937) and 
Maeda and Fujita (1989) on 6H-SiC who concluded that the BDT temperature o 
this material is above 800-1000'G The activation enthalpy Q (see ^bon (6) Owa 
determined from the slopes of the straight lines m figure 2 and found to be 
aWl 7 eV (in the T > T B regime). Using a value of n « 3, determined or, a <W- 
£ ttSctt * SamanTand Pirouz 0W ^«^ 
few strain-rate jump tests on the new 4H-SiC material by Demenet et al. (2000), gives 



OCT- 14-2002 11:42 



Linda Hall Library 



816 926 



1214 



P. Pirouz et al. 

115 9V __B__L(° C > 




14.0 



6.0E-04 



7.0E-04 8.0E-04 

i/r<m 



9.0E-04 



Figure 2. Plot of In (r y ) versus \/T for 4H-SiC deformed at strain rates of 3.6 x 10" 5 s" 1 (O) 

. and 2.6 x .10"° s ' (#). 

n ~ l 5 and 1 8eV at e « 3.6 x 10' 5 and 2.6 x lO'V 1 respectively, resulting in 

fa K X 0.3eV for the former strain rate Ai/ g * 5.1 ±0.3eV for the latter 

^Following the compression experiments, extensive TEM of the deformed samples 
was performed. In figure 3, the dislocation configurations in two of the samples 
deformed at a temperature higher than T c are shown: figures 3 (a) and (^correspond 
to deformation at T = 1200°C and i « 3.6 x 10" 5 s" 1 , whereas figures 3 (c) and (<*) 
correspond to deformation at T = 1280X and e « 2.6 x 10"< ^' 1 These are repre- 
sentative of the microstructure of many samples deformed at T > T c for both strain 
rates and show that deformation of 4H-S1C above the transition temperature gen- 
erates a large density of b-±<1120> dislocations dissociated into 5 (1010) leadmg- 
trailing partial pairs on the (0001) slip plane. It is also noteworthy that a Burgers 
vector analysis reveals that most of the dislocations have a screw character, presum- 
ably indicating that the mobility of non-screw dislocations is higher m this material, 
something that is consistent with other tetrahedrally coordinated materials. In con- 
trast with the samples deformed at T > T c , there is a very low density of dislocations 
generated by deformation at T < T c and these are predominantly 30° single leading 
partials all with the same Burgers vector b, = J (1010). Examples of these are shown 
in figures 4 and 5 for deformations at T = 1000 and 840°C at the higher 
(e«36x 10- 5 s" 1 ) and lower (e « 2.6 x 10~V l ) strain rates respectively. Figures 
4(a) and 5 (a) show the dislocation configurations on the (0001) slip planes; with the 
reflection g = 1120 used for imaging, the stacking faults behind the partials are out 
of contrast. On the other hand, in figures 4(A) and 5(b), the reflection g = 1101. is 
used to show the stacking-fault contrast. The appearance of overlapping stacking 
faults indicates that the leading partials are on parallel slip planes. In addition to the 
strain contrast experiments, the cores of a number of single partial dislocations, 
produced by deformation at temperatures below T e , were identified by the technique 
of LACBED. Examples of LACBED patterns for one of the dislocations in figures 
4(a) and 5 (a) are shown in figures 4(c) and 5 (c), respectively; in all cases examined, 
the core of the leading partial has turned out to be silicon rather than carbon. 
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Figure 4 (a) Bright-field micrograph of dislocations in 4H.SiCdefomied at T = 1000 C and ■ 
F ^iXlO- 5 s-' showing single leading partial ^locations; 0> same region as (^,- 
KJing ov«la P ping stacking faults; (c) a LACBED pattern of a single leadmg partial 
• in (a). . -\ 
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Figure 5 



• t»\ Rri^ht field microeraph of dislocations in 4H-SiC deformed at T = 840°C and 
~ ft SSi leading partial dislocations; (b) same region as (b 

l^JLtftf^A* (0 a 8 lWD pattern of a single leading P a rtl al 

in (a). . 



§3. Discussion 

Similar to figure 2, plots of the critical yield stress versus temperature for all the 
other five compound semiconductors mentioned previously (InSb, InP, traAs 
(Suzuld ,^// 1999b), GaP (Edagawa et al. 2000) and ^ Pj^ 
1998)) show an abrupt transition at a cntical temperature T It is possible that ^ 
transition exists in aU the tetrahedrally coordinated materials and signifies , a change 
in the deformation mode in the material. Also, intrigumgly, in every ™ 
tion temperature in the r Y (r) plot appears to coincide a PP^ a ^^ h ^ B k f ] t 
temperature r BDT of that material, taking due account of the strain f*. S ™^ 
al (1999a,b) explained the change in the deformation mechanism at T c m the com- 
pound semiconductors that they studied by a change in 'the slip plane (shuffle to 
glide) on which dislocations glide. Thus two regimes for dislocation motion in these 
materials were proposed: perfect screw dislocations moving in the. shuffle set at low 
temperatures versus dissociated dislocations moving in the glide set at high tempera- 
tures. These arguments were based on some compelling experimental evidence 
including the observation of wavy slip lines on the surfaces of samples defonned 
at low temperatures (indicating cross-slip, something that is only possible for perfect 
screw dislocations), small activation volumes at low temperatures, and a temperature 
dependence of r Y that is more consistent with kink pair formation in perfect screw 
dislocations rather than in partial dislocations. The observation of dissociated dis- 
locations by TEM in their low-temperature-deformed samples was attributed to 
post-mortem relaxation of. stress and pressure after deformation. In this way they 
assumed that, although perfect screw dislocations are responsible for the deforma- 
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u seriate into their equilibrium width after the 
tion at ^ t-P«;^^ (1999a,b) on the low 

deformation has stopped. The ugg • ^ huffle lane is consiSte nt with the 
temperature motion of per^^^ e , fl , (1994 ). Ren «r ^- (1995) and 

work of Kaxiras and Duesbery , m) calcu i ate d the Peierls Stress tor 

- Duesbery and Joos (1996). Thus Joo s et ££££ model with generated 

shuffle and glide dislocates ^ meta0 ds (Kaxiras and Duesbery 

stacking faults, determined ^^^Sgmtude higher for glide partial dis- 
1993), and found values nearly an order ■ ot magni by ReQ w 

ocations than for shuffle ^^^S^ the empirical Stillinger-Weber 
a/. (1995) employing atom f ^^^hemeen theoretical predictions of perfect 
( 995) interatomic potentials. ^^^oLcl^ly dissociated disloca- 
huffle dislocations and experimental <**g£™ afld Jo6s (1996) on the basis that, 
tioDSontte^pl^^^^^^^ieruJ. where there is significant 
at finite temperatures, in ^J^^^^^t^ir^^ ^^ 
dislocation-lattice coupling ^^^^^cuitfid the kink pair fonnation 
by ldnkpaknuoteitionaiidiiugra^n^tft ^ ^ ^ exCept 

enthalpy using a ^^7"^^^ enthalpy for a. (30-) partial i,s,g- 
very high stresses (above GPa^ ^ a ^a ^ indicating that slip on the 
nificantly lower than .that for a shuffle ^ s ^^ th experime nts. . 
glide planes is indeed preferred, in ^^f™^ ^describe this transition m the 

IxTthe following, we ^^^^ e t ai 1999a,b, 2000). 
yield stress and try to relate it to the BL> l (fir ^ 

3.1. 

As shown in the previous sectrot, our £EM results she of only 

of the SiC samples that were deformed ^below 1 o ^ ^ ^ ^ 
Si(g) leading partials whUe * e ™^ If ™ 

sisted predominantly of diss <«^^J^S deformation of such a tctrahed- 
n ow consider the product**! ^^^^Srway to produce dislocations in a 
rally coordinated crystal. In ^\^ S f Z & \ he Frank-Read mechanism, 
crystal that contains appropriate souses is y low tempera ture S 

^wever, if only P^,^-^, 6 ^^^), then such a source will 
(by the operation of the leading V^^^^^ once . Subsequently, 
not be a multiplication site for ^^ fr ™*^ ; he sarnp le or else by nucleating 
the stress will have to be relieved «^^ fr £S within die bulk), 
dislocations from the sample activation is often thought to 

In the energetics ^^ r ^^S to be that the available thermal 
play a negligible role. The ^/"^^Sd with the large enthalpy barrier 
energy kT at a temperature T is to^^a^ the taction rate theory, ther- 
Atf n required for dislocation no ^jnJt term* o ^ vibraie 

mally activated If^Sp^JS^^^ exp(-Atf n A^) of 
o vlrS ttS^li? A^dingly, the frequency „ with which such an 
activation will occur will be: 

, **r\ (7) 



„ = vc exp ^- J , 
, = (1 /„,) «p (Aft/Jtr). For djslocation nucleaoon, 
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i /mnzrr ur^Mer\ that the time i for a successful operation 
^&**"^ X ^*^l^$i for an activation enthalpy of 
SSJn St£S ^oc^d^e.uation (2)), the free energy for 
S^o n is also probably a functio n of the ^pW 
approximation, the variation in the free-energy change AG n <<* «^ J™£ 
AH ) with the applied shear stress r may be assumed to be linear , that is the effective 
fnthalpy barrier Atff may be assumed to be reduced by the applied stress r accord, 
ing to the following equation 



A/£ ff « A# n - arK*. 



(8) 



In this equation, ar is the stress concentration at a heterogeneity (e.g. a scratch at the 
uiace and" is the activation volume (- -(9H,/dr) r ). Thui. In the presence of a 
Sear stress, the enthalpy barrier AH n for dislocation nucleation 
amount arK* to AH* . Accordingly, the probability of overcoming the barrier 
becomes exp \—{AH B — arV*) /kT] and the frequency of successful dislocation 
SI Events will be „ - i'exp - «r7*)/*I1. Hence, the mean time t 

during which a successful activation occurs will be 

(9) 



1 fA H n -arV*\ 

=ir xp l — kr — y 



Based on the experimental results, we shall assume that, in the above equation the 
leading and trailing partials have different activation enthalpies for dislocation 
nucleation, AJ* B and AH x a respectively, and that Atf]> < A/4 In addition, in equa- 
tions (8) and (9), the resolved shear stresses on the leading and trailing partial 
dislocations will be different because they have different Burgers vectors. 

Since our experiments were performed on hexagonal SiC> we shall use this mate- 
rial as an example, but the arguments are essentially general and can be applied to 
other semiconductors. Figure 6(a) shows a schematic diagram of the orientation of 
the parallelepiped sample used in our compression experiments. The sample is 
oriented for single glide, that is an orientation in which the (0001)(21 10) slip system 
is primarily activated. The dislocations most probably nucleate heterogeneously as 
half-loops from the crystal surface; an example of one (denoted as A) on a particular 
(0001) slip plane is shown in figure 6(a). For comparison, the nucleation of a 
dislocation (full-)loop (denoted as A') on the same (0001) slip plane within the 
sample is also shown. In fact, nucleation of four types of loop on a (0001) sUp 
plane may be envisaged, as illustrated in figure 6(b). In cases A and A , only the 
leading partial has nucleated while, in cases B and B', nucleation of the leading 
partial has been followed by that of the trailing partial. 

The formation of a dislocation loop presumably takes place by, first, the nuclea- 
tion of the leading partial, followed by the nucleation of the trailing partial. 
Assuming that intrinsic faults are favoured over extrinsic faults, then the order of 
nucleation of the two partials can be determined from the Thomson tetrahedron 
(Hirtb and Lothe 1968). It is important to note that once the leading partial forms, it 
creates a stacking fault behind it that changes the state of the lattice for nucleation of 
the trailing partial. The stacking fault of course adds an extra stress term j/b to the 
external applied stress that should make it simpler for the trailing partial to nucleate. 

It was mentioned above that the assumption has. been made that Ai/ n and Atf fl 
are different and that AH l n < A^. In fact, there is some independent evidence for 
this hypothesis in GaSb and GaAs (Ning et al. 1995) and also in silicon (Ning and 
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Huvey 1996). There is, however even - , ^£5££* 
leading and trailing partiais m example Wessel and 

their mobilities (for the case of silicon and PJ«» ^ this difference 

Alexander (1977) and Alexander e t ai ^^L^portan.* is the mobility 
ceases as the J^SSSSd screw dislocation where it 

difference between the two 30 partiais 01 a - u and ger manium, the 

is found that, even in ^f^^^S^^ P-**- ™ S ^ 
fading 30° partial *f * are also 

■ ihat the activation ^J™^*^* s .nobility and nucleation differ- 
different and that AN 1 , < Atf g; probabry have a similar origin, 

eiices of the two partiais are not really ckar, ^ *J £ r ° ^ ^ formatloa of the 
One obvious difference » the. presence o a suclan ^ nUC leation of the 

leading partial that provides a l^^S of a dislocation in a 

trailing partial. However, it * ^^f^' 1 ^^^! nucleation in a perfect 
faulted region of the crystal wdl be J*^^™^ ^ stress m that 
crystal. Moreover, the presence of a stackin ^ P her sibiUty is that 

should actually help the nucleation of the tra lrn ^ arU * ^ £ ealiorL of p0 int 
the fonnation.(and glide) of the leading 5 partml may ^ ^ 

defects in its wake (i.e. in the faulted parti al is nucleated, 

aificant effect on the local environmen 1*^™^?^^**^ 
These point defects need not £ course be "J^^^ interaction of point 
antiphase defects (Hirsch 1980, Heggie and J ™ ^> back to the early 

defeats with dislocations and it, f*™*?^™™^ matnals with a high 
work of CeUi et al. (1963) and ^^^^^viet scientists perfonned 
Peierls barrier. Also, m the late 1970s and 198US tne ao 

many experiments on silicon and g = £S^a,, thus 
swept by a moving dislocation is ^^^Z f xample fiondarenko « «/. 
affecting the motion ot sub se ^ e f^^ 

1981). Some important ^^Z^^^%o^aa of point defects by 
of Kisielowski-Kemmench ( 99^ tto JT^^ of locations in semicon- 
screw dislocations m silicon. In fact, m hi s cla^ ^ dis i 0 cation-point- 
ductors, Alexander (1986) emphasizes ^^^^^ of dislocations in 
defect interactions is essential to m t rad !^^^ t U^ usually annihilate 
these materials. It is also important to realize that the point oeiec 
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(or anneal out) at higher temperatures and their interaction with dislocations 

becomes less effective. , . r *\.~ \^^a 

In the case of compound semiconductors, where the atomistic cores of die lead- 
ing and trailing partials may be different, the core nature is another factor that miut 
be taken into account in considering dislocation nucleation. Thus, m the case , o S C 
the partial with the silicon core may have a smaller activation enthalpy than the 
partial with the carbon core. As was mentioned before, half of a partial loop (with 
any Burgers vector) has a silicon-core while the other half has a carbon core (figure 
n We suggest that, in the case of SiC, the Si(g) half-loop of the leading partial 
(whose Burgers vector is determined from the Thomson tetrahedron) is preferentially 
nucleated compared with a completely C(g) or a mixed Si(gKXg) half-loop (figure 

^Assuming that, at T < T c , only the leading Si(g) partial (Burgers vector b,) is 
nucleated, then the sequence of events is illustrated in figure 7. After its nucleation 
from the sample surface (figure 1(a)), the half-loop expands and, because of the 
different mobilities of the segments, the semihexagonal shape of the half-loop is 
distorted. Assuming that the 90° Si(g) partial has a higher mobility than the two 
adjoining 30° partials, the configuration of The half loop at a later stage becomes, 
something like that shown in figure 7 (b). In figure 7 (c), the 90" partial has moved out 
of the specimen and the lower 30° segment (which has been assumed to be more 
mobile than the upper 30° segment) is in the process of moving out. The final out- 
come will be a preponderance of 30° Si(g) dislocations left on different (and parallel) 
(0001) planes. Note that only a few partials can exist on each particular (0001) plane 
(figure 6(b)) because each partial drags a stacking fault behind it and the initial 
source (AD in figure 7 (a)) is deactivated as a dislocation source. 

Now consider the nucleation of a perfect dislocation half-loop from the sample 
surface (B in figure 6 (b)) at T > T c . The atomistic detail of the half-loop B in figure 
6(b) is shown in figure 8 where the cores of all the half-loop segments are shown 
reconstructed. The sequence of events after nucleation, that is the expansion of the 
perfect dislocation half-loop on the (0001) slip plane, is illustrated in figure 9. Here, a 
dissociated dislocation half-loop, ABCD, nucleates from the surface (figure 9(a)) 
and expands (figure 9(b)) until the mobile segments run out of the sample (figure 
9 (c)), leaving only the slow dissociated screw dislocation on the (0001) slip planes. In 
this figure, additional assumptions have been made that mobility of the upper 
90 6 Si(g)-30°Si(g) segment is higher than that of the lower 30°Si(g)-90 o Si(g) segment, 
which in turn is higher than the topmost 30°Si(gH0°C(g) sctew segment. The final 
configuration will be a preponderance of 30°Si(g>-30 o C(g) screw dislocations parallel 
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Figure 7. Nucleation and expansion of a leading partial dislocation half-loop from the 

surface of the crystal at T < T e , 
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to the [2110] direction in the deformed crystal. Clearly a better knowledge of the partial 
^nb« tie i essential to draw clearer conclusions. Multiplication of dislocations from, 

« Jto totty of (perfect) dislocations in the crystal exceeds about 10 -10 cm . 
"rtS^miSti^tici-ntfhei than nucleation, will take over as the dominant 
5332 S basing the dislocation density for deformations above r e . 

3 2 Relation between T c and Tbdt 

Tie next problem to consider is the possible relation between the transition in the 
T (T) tdoT aSd the BDT temperature. When a load is applied to a crystal, it » 
3 2l U Sfottrt elasticaUy whereby elastic energy is stored in the crystal This 
S Energy is released if the load is taken off. If, however load.ng ,s 
continued at a ratf r, and the stored energy continuously increases here are two 
ma^nTompeting ways for the crystal to relax and decrease its overall elastic energy. 
S^e way^is by brittle fracture of the crystal whereby the energy is consumed in 
creating two new surfaces as a crack forms and propagates. 

Se second way for the release of the stored elastic energy is to move pre-extstmg, 
dislocations or to nucleate and move fresh (as well as pre-existing) dislocation^ The 
JcurrTnee of brittle fracture versus plastic yielding of a crystal is determined by the 
competition between these two ways of releasing the stored elastic energy. In the case 
of crack propagation, the Griffith criterion can be written in the form 



(10) 



where v is Poisson's ratio, E is the elastic modulus, 7s is the surface energy and K eir is 
the effective stress intensity factor given by (Thomson 1986) 



Kffr == -Kapp &i 



with 



appi 



(11 




Figure 8 Schematic illustration of the atomistic detail of dissociated (leading jailing) 
^ dislocation half-loop nucleated from the surface of the crystal at T s 7 C . 
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Figure 9. Nucleation and expansion of a dissected (leadin^traiUng) location half-loop 
s from the surface of the crystal at T > T z . 

~ hrine the annlied tensile stress and a the crack half-length. & in equation (11) is 
&J^dfSJS^Ud«d by the compressive stress of the dislocaUons 
afou^dTthis results in a reduction of the applied tensile ™ *'^ tal 

In the above conventional picture, brittleness 
is considered as the competition between shearing of ^^ty**^ 
motion and rupture of the bonds at the tip of a microcrack ™ to the crys taL" The 
general vrew is that the stress needed for plastic yielding (i.e. the y-l\ stress r Y ) 
lonely temperature dependent and decreases rapidly with increasing temperature 
w^e (tensile) stress^ 

shear stress r F ) has a much weaker temperature dependence (Kelly et al 1967). As a 
esuk r'^) and r Y (T) intersect at a critical temperature r BDT , usually known a 
the !bDT temperat ure, where r F < r Y at T < T BDT (the brittle regime) and r Y < r F at 

T Paper (as illustrated in figure 10) the B^at least 

in tetrahedrally coordinated materials, is determined by the coro P^ mo \^^ n n ^ 
nucleation and" propagation of leading partial dislocations versus the nuefca, onand 
propagation of perfect (total) dislocations (i.e. that of the trailing pa rtia 0 A» sh«£ 
by equation (9) the nucleation. and propagation of dislocations (part lal or ^perfect) 
are temperature and stress dependent. Thus, the rupture ot atomic bonds arte* the 
picture only indirectly. It is argued that, in a constant-stram-rate «*Pf™^ ™ 
temperature and stress are not sufficient to nucleate the trailing partial (and thus the 




Figure 10. Schematic illustration of dislocation nucleation at a crack tip under mode 1 

loadingt (a) at T < T B ; (/>) at T > T c . 
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„a .tress will eventually reach a sufficient value 
perfect) dislocations, then the increas* propaga tion. In such a regime, 
? F t o rupture the bonds at the crack ^*£ at th / crack tip (as well as id other 
£ding partial dislocations may well jw^ate^t ^ . & ^ by 

^oufable sites) but, if at ^.P"^^^ shut off and stop operating 
nucleation of the trailing partial *^ J*^^ soon stop and plastic dcforma- 
The consequence is that location nuckat^n dislocationS gene rated before 

tion is limited to the strain P rodu ^ w Me the work done on the crystal by the 
he sources stopped operating. ^^^^ energy until the critical strain 

external load continues ^T^J^U^ * ™ d fractUreS " * J?; 
energy release rate is ^/ c ^p^ temperature dependence ofay (or r*) ate 
details of bond rupture at the cracK up o 

"relevant to the problem of the BDT. ^ . g ^ temperat ur C depen- 

t^aS £ S; total) ^^J^^^S. r (rela,d - ; he 
the intersection of t> (7) and T« (T) f ^-^^ factor) is applied at a rate ot A 
Sear stress r by K /TV, ^^VZc^ul (dislocation) nucleation event wil 
the question to be asked is whether a successtuH ^ ^ ^ th 

occur during the time period m which K * n) . Note that, as A' (i.e. r) 

*S LJ* ^ctor X lc at which the cry t* fob (fig- ^ nudeation . decr ea S es. 

ncreases, the effective c"uo^ 

Thus the BDT is determined by a c°mP ei or _ a /jJAP 

Vh ± K reaches K lc first and catastrophic tauu dis locanon. 



Thu7 th; BDT is determined by a or - o/«» 

events- either K reaches K lc first and catas ^^ (tnriUng partial) dislocanon. 
«s%cienti y to result m the | nude£ the condign for 

The time taken for K to reach X te is (*.c *oJ/ 
BDT is that 

, *- jr it 7W to nucleate dislocations. 
Let the stress intensity factor reach a value K = Ay at to 

Thus r = (*r - *>)/*• ^ Vine 




'time 



Figure 



11. Schematic plot of K versus time ma 
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(13) 



where I = °PV* « a constant. Tims, a plot of ln(*) versus l/T„ 
wnere o cv iJCvVk and an intercept In f 0 (*.y - KoJI- 5m< * 

straight hue with a slope (AH. <^""* ._ ^ a slope equa l , 0 that of 
experiments show that a plot of In (K) versus 1//bdi na ^ ^ 1?87 

^^^l£^f« nuclea.™ of — 

given by A*f =^"^ — a ^^J**£X one 

S"*-:" aXSon enthalpy front the vafcesofn a t *e 

obtains 

. In f 3.6 x 10- 5 /2.6 x 1Q- 6 ) _ t 
(1/1373 - 1/1303) 

^rnenS aS^nei to obtain the shift of T a, ^ ° £ 

, i i < ^ iir\ vAfQii^ 1/T is indeed linear, if so, a better vaiue ui mo 

"T*H t . Shjar experiments on other semiconductors wonld also be hrghly desir- 

able. 

3 3 Incidence of perfect (or narrowly dissociated) screw dislocations during 
very-loy> temperature (T « T c ) deformations 

A number of experiments have .shown that when high-quality (very-low-dis toca. 
tion^^Jm^faicl^ crystals are deformed under hydrostauc pressure there 
fs a range of temperatures T BDT > T > r„. 

r-Saizuki ^ al 1999a b Edagawa et at. 2000). This is so much so that, if the crystals 
b teTot« 

small pieces during decompression of the sample cell. An example is GaAs that was 
2£d££ — y britL between 300 and 400 K even when 
sure was increased up to 1.2 GPa; consequently, it has been ^j^**?*^ 
thin foils and to perform TEM experiments on the samples M°™?. m J^™% 
(Suzuki et al. 1999a). However, interestingly, below this range that is a ^mpera^ 
ures less than T u , it again becomes possible to deform the c ^^^^^ 
stresses and under hydrostatic pressure. Thus, Suzuto W d. (1999a ) ^d^form 
GaAs at 300 K under a confining pressure of 400 MPa! TEM of samples deformed at 
these very low temperatures show a microstructure consisting of P^onunantly 
perfect screw dislocations. At first sight, this appears to be contrary to the argument 
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, tw nnlv leading partial dislocations may be nucleated at 
put forward in this paper that ive activation enthalpy for 

T < Tbdt- However, according to . equafi ton _ ^ fey ^ appUed 

dislocation nucleation given by A* ^ ^ (ovgr Q 5Gp } 

resolved shear stress r. Thus, it may d j sufficiently to make 

applied at temperatures below fu.low toa^a^ resultrng in leading-trailing 
Ration of the trailing partial fifi^^ £ also supported by 
pairs of partials, that is perfect -^^^^^^ thai the n«deacion 
theoretical calculations of Br0c ^ ard ef a j ( a \ 9 ;Xiently high applied stresses, both 
of dislocations is stress at moderate.y 

partials can nutate f - a ^XSons are nucleated sequential!, .(leading 
low temperatures (T < T a ), partial ui temperatures (T Ju). 

foUowed by trailing) on the ^fP^^^Ls becomes very high, it 
where the shear stress required o nudea d ^m*™ to nucleate 
becomes energetically favourable for ^\^°^ s bery. the screw nature of 
and propagate on ^J^^ZS^ «f non-screw segments 
the dislocations ,s probably due to ™ idl m0ve out 0 f the crystal and 

of a half-loop which, under the very high stiessev rap y 
leave the much slower screw segments m the sample (as m ngure ; 

84 Conclusion 

New deformation experiments on ^^X^ ^^' 
temperature in the ^j^?^^^ ^ " ^ 

Similar to the lowering of ^.^J^^l also lowered with decreasing strain 
variation in the yield stress with tempe ature . also ^ ^ 

rates. The decrease in T c is about 70 C trom 1100 C at ^ ^ 
at , = 2.6 x 10- • TEM experimen^ve ^^J^^ single 
locations produced in the samples ^^J^T^L of S i(g)-C(g) dissociated 
leading partials with a ^^^^^^^JSt^A model has 
screw ^^^J^^ZSSZ^^ T C with a change in the 
been presented that asso «ates tne u dislocations are nucleated 

mod e of deformation. In this mode , ordy lading _P ^ wtent 

at r < r ei which glide on parallel slip planes but ^tribute y ^ 

to the stoning of the crystal, espeaaUy sin ce ^Xltion'on the other 
become inoperative after they emi oriy one ^ from sourceS iu the 

hand, at T > T c , full dislocations ^ n * n »^ of the 

crystal and their glide ^^^^^^^Jt c in the yielding 
crystal. Additionally, the model relate s * at P ^ 

of the crystal to ~ ^ *t very low temperature, iu 
occurs. The formation of rfeCt h lied stre sses that lower the activa- 
semiconductors is attributed toto ^'g dislocations. Deformation 
tion enthalpy both at moderately 
experiments m the brittle regime oi vano characterization of the result- 
low and very low temperatures, light on the mode 
ing dislocation "f of theTe materials, 
of dislocation nucleation and the transition m 
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